Spectral spacing correlations of nucleus-nucleus collisions at high energies 
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We propose a novel approach to the analysis of experimental data obtained in relativistic nucleus- 
nucleus collisions which borrows from methods developed within the context of Random Matrix 
Theory. It is applied to the detection of correlations in momentum distributions of emitted particles. 
We find good agreement between the results obtained in this way and a standard analysis based on 
the two-pair correlation function often used in high energy physics. The method introduced here is 
free from unwanted background contributions. 
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There is currently an enormous effort underway to 
detect signals of possible transitions between different 
phases of a composite system produced in high energy 
nucleus-nucleus collisions (cf [HQ). It is anticipated that 
in central collisions, at energies that are and will be soon 
available at SPS(CERN), RHIC(BNL) and LHC(CERN), 
the nuclear density may exceed the density of stable nu- 
clei by an order of magnitude. Under such extreme con- 
ditions, according to a generally held beliefs, the final 
product of a heavy ion collision would be a composite 
system that consists of free nucleons, quarks and a quark- 
gluon plasma. Various methods have been proposed to 
identify possible manifestations of such a quark-gluon 
plasma. Often though, results based on such methods 
are sensitive to assumptions made concerning the back- 
ground measurements and mechanisms included in the 
corresponding model. In addition, the identification of 
the quark-gluon plasma is made more difficult due to a 
large multiplicity of secondary particles created at these 
collisions. The formulation of a reliable criterion for the 
selection of meaningful signals is, indeed, an important 
objective in relativistic heavy ion collisions physics. 

In a preliminary report [3|, we suggested that tools 
from Random Matrix Theory (RMT) [4j] might be useful 
in illuminating the presence of correlations in the spec- 
tral (momentum) distribution of secondary particles pro- 
duced in nucleus- nucleus collisions at high energy. The 
RMT approach does not depend on the background of 
measurements and relies only on fundamental symme- 
tries preserved in heavy- ion collisions. Furthermore, the 
larger the multiplicity, the better the applicability of the 
RMT tools, and thus, its predictive power. In the present 
paper, we demonstrate that the RMT analysis is very 
sensitive to spectral spacing correlations present in the 
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nucleus-nucleus collision data, more so than the standard 
tools used for such an analysis. 

Here, we make use of the experimental data that have 
been obtained with the 2-m propane bubble chamber of 
LHE, JINR [5], [||. The chamber, placed in a magnetic 
field of 1.5 T, was exposed to beams of light relativis- 
tic nuclei at the Dubna Synchrophasotron. Nearly all 
secondary particles, emitted at a Air total solid angle, 
were detected in the chamber. All negative particles, ex- 
cept those identified as electrons, are considered as ir~- 
mesons. The contamination from misidentified electrons 
and negative strange particles does not exceed 5% and 
1%, respectively. The average minimum momentum for 
pion registration is about 70 MeV/c. The protons were 
selected by a statistical method applied to all positive 
particles with a momentum of \p\ > 500 MeV/c (we iden- 
tified slow protons with \p\ < 700 MeV/c by ionization in 
the chamber). In this experiment, there are 37792 12 CC 
interaction events at a momentum of 4.2A GeV/c (for 
greater discussion of the details see @) containing 7740 
events with more than ten tracks of charged particles. 

The basis of our approach derives from RMT Q , which 
was originally introduced to explain the statistical fluctu- 
ations of neutron resonances in compound nuclei 0] (see 
also Ref. Q). There the precise heavy-compound-nuclear 
Hamiltonian is unknown or rather poorly known, and 
there is a large number of strongly interacting degrees of 
freedom. Wigner first suggested replacing it by an en- 
semble of Hamiltonians which describe all possible inter- 
actions @ . The theory assumes that the Hamiltonian be- 
longs to an ensemble of random matrices that are consis- 
tent with the fundamental symmetries of the system. In 
particular, since the nuclear interaction preserves time- 
reversal symmetry, the relevant ensemble is the Gaussian 
Orthogonal Ensemble (GOE). Whereas if time-reversal 
symmetry were broken, the Gaussian Unitary Ensemble 
(GUE) would be the relevant ensemble. The GOE and 
GUE correspond to ensembles of real symmetric matrices 
and of complex Hermitian matrices, respectively. Besides 



the general symmetry considerations, no other property 
of the system under consideration is taken into account. 

If one supposes that the momenta distributions of sec- 
ondary particles produced in nucleus-nucleus collisions 
may be treated in analogy with eigenstates of a com- 
posite system, just like the eigenstates of the compound 
nucleus, then the statistical analysis methods introduced 
by Dyson and Mehta can be applied to the LHE collision 
data [7j. The difference between energy and momentum 
is not essential for pions, and we assume that the proton 
mass does not significantly affect the correlation function. 
Note also, that here we are dealing with the momentum 
distribution in the target rest frame only, postponing its 
comparison to that in the center of mass frame, which 
is more natural for description of interaction. Based on 
this supposition, the ordered sequence of "energy levels" 
{Ei}, i = 1, N comes from the momentum distribution 
and it has an average density of states denoted p av (E) . 
From this sequence a new one is obtained by the unfold- 
ing procedure of the original spectrum {Ei} through the 
mapping E — ► x 
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FIG. 1: (Color online) Cumulative spacing momentum distri- 
butions for different regions of measured momenta: a) 0.1 < 
\p\ < 1.14 GeV/c (solid line); b)1.14 < \p\ < 4.0GeV/c 
(dashed line); c)4.0 < |p| < 7.5 GeV/c (dot-dashed line). 
The solid line exhibits the Poisson distribution, while the dot- 
dashed line is close to the Wigner-Dyson type distribution. 



xt= / P av{E')dE' = / dx' = {(Ei), i = l,...,N 
Jo Jo 

(1) 

Here, ((E) is the smooth function giving the mean num- 
ber of eigenvalues less than or equal to E of the ex- 
act eigenvalue distribution N(E), which is often referred 
to as the staircase function due to its appearance (see 
Refi). The smooth part ((E) can be determined ei- 
ther from semiclassical arguments or by using a polyno- 
mial/spline interpolation for the exact staircase function. 
The momentum distribution data (see Fig.l in Ref.3) 
were approximated by a polynomial function of sixth or- 
der and the distributions of various spacings Si from the 
7740 events satisfiy the condition of x 2 per degree of free- 
dom is less than unity. 

The effect of the mapping is that the sequence {x{\ 
has on average a constant unit mean spacing (a con- 
stant unit density), irrespective of the particular form 
of the function ((E) The spacings are defined as 

Si = Xi + i — Xi between two adjacent points and they 
are collected in a histogram, which gives the probability 
density. If the "events" {xi} are independent, then the 
form of the histogram must follow p(s) — exp(—s) known 
as the Poisson density. On the other hand, if the lev- 
els are repelled (anticorrelated) as in the GOE, the den- 
sity is approximately given by the Wigner surmise form 
p(s) = ^sexp(— f s 2 ). Interestingly, the spacing proba- 
bility density approximately follows the Wigner surmise 
for high energies, whereas at relatively low energies the 
corresponding spacing density is maximum at the origin 
and nearly the Poisson density [3j . In an eigenspectrum, 
the Poisson density arises where there is a dominance of 
many crossings between different eigenenergies, whereas 
the Wigner surmise reflects the tendency to avoid cluster- 
ing of levels. In turn, the crossings are usually observed 
where there is no mixing between states that are char- 



acterized by different good quantum numbers, and the 
anticrossings signal a strong mixing due to a perturba- 
tion brought about by either external or internal sources 
(cfRefsHi). 

The transition from one probability density to the 
other has been used, for example, in nuclear structure 
physics to study the stabilization of octupole deformed 
shapes and transition from the chaotic to the regular pat- 
tern in the classical limit 10]. Therefore, such an analysis 
can provide the first hint of some structural changes at 
different parameters of the system under consideration, 
in particular, in different energy (momentum) intervals. 
Figure Q] shows the integrated momentum spacing den- 
sity for experimental data, demonstrating a gradual tran- 
sition from a Poisson-like density toward a Wigner-like 
density with the increase of the absolute value of the mo- 
mentum distributions. 

In order to elucidate the degree of correlations for a 
stationary spectrum with unit average spacing Dyson in- 
troduced the k-level correlation functions 

R k (x\, ...,Xk) = (jy-fc)! / p k(xi, ...,x N )dx k+ i...dx N 
1 < k < N, (2) 

where Pk(x\, XN)dx±...dxN gives the probability of 
having one eigenvalue at x±, another at X2---, another 
at x 7v each within the interval {xi,Xi + dxi}. By inte- 
grating Pk(x±, xn) over all variables but one, in the 
limit N — > oo, one obtains the ensemble averaged density 

p(xi) = J ... J Pk(xi,...,x N )dx 2 ...dx N (3) 

which is normalised to unity. From Eq.@ it follows 
that i?i(xi) = Np(x\) and Rk(x\, Xk)dx\....dxk is the 
probability, irrespective of the index, of finding one level 
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within of each of the intervals [xi,Xi + ds]. From the 
above definition it follows that R\{x) = 1. With the aid 
of the definition J2|), by integrating R k +i one obtains 



R k+ i(xi, ...,x k +i)dx k+ i = (N-k)R k (xi,...,x k ) (4) 
It is difficult to work directly with the R k functions. One 




FIG. 2: (Color online) Number variance S 2 (I/) for three dif- 
ferent region of the experimental spacing momentum distri- 
butions: a)0.1 < |p| < 1.14 GeV/c (solid line); b)1.14 < |p| < 
4.0GeV/c (dashed line); c)4.0 < \p\ < 7.5 GeV/c (dot-dashed 
line). 

important and more convenient measure of correlation 
that was introduced is based on the number statistic n{L) 
which is defined to be the number of levels in an energy 
interval of length L. Since the spectrum was unfolded, 
the average number statistic (n(L)) = L is independent 
of the spectrum. However, the variance of n(L) 



Z 2 (L) = ([n(L)-(n(L))} 2 ) 



(5) 



does depend on the spectrum considered. For the Poisson 
density (see [H) 

Y?{L) = L , (6) 
and for the GOE, the exact asymptotic expression is 



7T^ 



ln(2nL) + 7 + 1 + ^[Si{nL)} 2 - ^Si(irL) 



- cos(2ttL) - Ci{2irL) + it 2 L 



' : Si{2nL) 



(7) 



Here 7 is the Euler constant and Si, Ci are the sine 
and cosine integrals, respectively. The number variance 
S 2 (L) calculated using the optimal implemention of the 
definition in Eq.© is shown in Fig.O Indeed, this quan- 
tity manifests the Poisson statistics (Eq. ([6])) for experi- 
mental spectra with a low momenta distribution. On the 
other hand, one again observes a clear indication on the 
presence of correlations for large momenta. 



For the analysis of fluctuations, it is more convenient 
to use pure k-point functions [12J 



Rk{L) 



R k (xi,...,x k )dxi...dxk (8) 



The function R k (L)/k\ gives the probability that an in- 
terval of length L (for small L) contains k levels. In RMT 
most emphasis has been put on the two-point correlation 
function i?2(a;i,^2) or density-density correlation func- 
tion. The two-point correlation function is the probabil- 
ity density to find two eigenvalues Xi and Xj at two given 
energies, irrespective of the position of all other eigenval- 
ues. The function R%{x\, X2) depends only on the relative 
variables s = x\ — X2- The variance S 2 is connected to 
i?2 through the following relation 



R 2 (L) 



i?2(xi, x-i)dx\dx2 



Z 2 (L) + L(L-l) 



2 / (L-s)R 2 (s)ds 
'0 

(9) 



The two-level correlation function R2 (xi,X2) determines 
the basic fluctuation measures related to Wigner's level 
repulsion and the Dyson-Mehta long-range order, i.e., 
large correlations between distant levels. Bohigas et al 
(1985) 12j provided a thorough analysis of level repulsion 
and long-range correlations (rigidity) for different corre- 
lation functions. To understand the distinct role played 
by level repulsion and long-range order in the momentum 
density, we compare our numerical results with analyt- 
ical expressions from Table 1 of ReflU for the Poisson 
ensembles (there is neither level repulsion nor long-range 
order) and for the GOE case (this ensemble exhibits both 
features). The two-point correlation function it^t^ii ^2) 
calculated from Eq.© is shown in Fig[3l Even though 
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FIG. 3: (Color online) Two-point correlation function i?2(s = 
L) for different regions of measured momenta: a)0.1 < \p\ < 
1.14 GeV/c (solid line); b)1.14 < \p\ < 4.0GeV/c (dashed 
line); c)4.0 < |p| < 7.5 GeV/c (dot-dashed line). The solid 
lines, denoted as P and GOE, display the characteristic limits 
for Poisson and GOE ensembles, respectively. 
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there are small deviations from Poisson (R 2 (L) = L 2 ) 
and GOE (R 2 (L) — 7r 2 L 3 /18) predictions, the experi- 
mental results for the momentum distributions reproduce 
surprisingly well both limits. 

The validity of the RMT analysis is confirmed by an 
independent analysis of the data with the aid of the stan- 
dard pair-correlation function (see, for example, Reffl3l 
and references therein): 



R{vx,V2) = <y 



d 2 a/dy 1 dy 2 
(da/dyi)(da/dy 2 ) 



1 



(10) 



Here, the quantity a is the cross section of the inclusive 

is the rapidity, which depends 



reaction and y = gfa ^p 1 



on the particle energy E and its longitudinal momen- 
tum P| | . The rapidity is one of the main characteristics 
widely used in relativistic nuclear physics (see [3, EH). 
The pair-correlation function manifests the difference be- 
ll 
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FIG. 4: Integrated two-pair correlation functions for particles 
obtained in CC-interactions (see the text) for different regions 
of measured momenta: a)0.1 < \p\ < 1.14 GeV/c; b)1.14 < 
\p\ < 4.0GeV/c; c)4.0 < \p\ < 7.5 GeV/c. 

tween probability density of two-particle events and the 



product of the probability densities of independent par- 
ticle events. It vanishes if the particle rapidities are inde- 
pendent. Figure |4] demonstrates the results for particles 
obtained in CC-interactions. For the sake of illustration, 
we integrate the function R{y\ 1 y 2 ) over one of the vari- 
ables, say yx, and consider the dependence on y 2 . For dif- 
ferent momentum distributions, there are three intervals 
of integration for the variable y±: a) for 0.1 < \p\ < 1.14 
GeV/c the function R(yi,y 2 ) is integrated in the inter- 
val -0.9 < j/i < 2.5; b)for 1.14 < \p\ < 4.0GeV/c it 
is integrated in the interval 0.5 < y\ < 2.4; and c)for 
4.0 < \p\ < 7.5 GeV/c it is integrated in the interval 
2.5 < yi < 3.5. There are, of course, some experimental 
errors which are not seen in the behaviour of the function 
R = f R(yi,y 2 )dyi on Figs|lk.,b due to the large data 
set. These errors do not spoil, however, the results in the 
region c), where experimental data on the multiparticle 
production exhibits an indication of the existence of cor- 
relations between particles in the region 4.0 < |p| < 7.5 
GeV/c. The physical origin of the correlations is beyond 
the scope of the present paper and will be discussed else- 
where. One observes that the predictions based on the 
standard pair-correlation function R are consistent with 
the predictions based on the RMT analysis. However, 
the RMT two-point correlation function magnifies the 
presence of correlations manifested in the standard pair- 
correlation function (compare Figsl3]and|4j:). 

Summarizing, we propose an analysis of relativistic nu- 
clear collision data based on ideas from RMT. The ap- 
proach is free from various assumptions concerning the 
background of the measurements and it provides reliable 
information about correlations induced by external or in- 
ternal perturbations. All these features make our pro- 
posal a quite promising avenue for the future analyses of 
data produced in heavy ion collision experiments. 
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